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Abstract 

We estimate the twist-4 spin-1 nucleon matrix element /2 in an instanton-based 
description of the QCD vacuum. In addition to the flavor-nonsinglet we compute 
also the flavor-singlet matrix element, which appears in next-to-leading order of 
the l/iV c -expansion. The corresponding twist-3 spin-2 matrix elements d<i are sup- 
pressed in the packing fraction of the instanton medium, p/R <C 1. We use our 
results to estimate the leading 1/Q 2 power corrections to the first moment of the 
proton and neutron spin structure functions G%, as well as the intrinsic charm con- 
tribution to the nucleon spin. 



1 E-mail: 

2 E-mail: 

3 E-mail: 



lee@tp2.ruhr-uni-bochum.de 

gockc@tp2.ruhr-uni-bochum.de 

weiss@tp2.ruhr-uni-bochum.de 



1 Introduction 



The structure of the nucleon as measured in deep-inelastic scattering is described by matrix 
elements of QCD operators of a certain twist (dimension minus spin). The moments of the 
structure functions at asymptotically large Q 2 are given by matrix elements of operators 
of twist 2, which have a simple interpretation as number densities of the quarks and 
antiquarks in the nucleon in the infinite-momentum frame. Power (1/Q 2 -) corrections 
to the asymptotic result are determined by matrix elements of operators of twist 3 and 
4. The latter describe either the correlation of the quark fields with the non-perturbative 
gluon field in the target, or quark-quark correlations. In the polarized case the twist-3 
and 4 matrix elements of lowest spin are that of the twist-4 spin-1 operator 

(p,X\iP na F^ f \p,X) = 2M 2 N f 2jS P, (1) 

and of the twist-3 spin-2 operator 

(p, \\ipf (y a F^ + 7 /3 F° 7 ) ip f \p, A) - traces 

2pVV - P>V - P y/ + (a++P)- traces! . (2) 
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Here ipf,ipf denote the quark fields of flavor /, 

= )f vp °F pu (3) 

is the dual of the gauge field strength (e 0123 = 1), p the nucleon four-momentum (p 2 = 
Mjv), A the helicity, and s = s(A) the polarization vector of the nucleon state, which 
satisfies 

s-p = 0, s 2 = -M 2 N . (4) 

The nucleon states are normalized according to (p, X\p', A') = 2p°(2n) 3 5^(p — p')5\\i. The 
matrix elements (0) and (|2|) can be taken either in a proton or neutron state. 

The matrix elements ([I]) and (|2|) together determine the leading power corrections to 
the first moment of the spin structure function G\ (the Bjorken and Ellis- Jaffe sum rules) 
|T|, |3j . The twist-3 matrix element d 2 also appears in the QCD expression for the third 
moment of the second spin structure function G 2 , where it is not power-suppressed relative 
to the twist-2 part and can therefore be measured with good accuracy [|J. By studying 
these matrix elements in theoretical models of the nucleon one may hope to get some 
insight into the transition from the region of asymptotically large Q 2 (Q 2 ^> lGeV 2 ), 
where the Q 2 -dependence of the structure functions is described by perturbative QCD, to 
the resonance region (Q 2 « 1 GeV 2 ) §. Eventually, this may help to understand also the 
discrepancy of the moments of the structure functions measured at Q 2 > 1 GeV 2 with the 
GDH sum rule for photoabsorption (Q 2 = 0) 0. 

The twist-4 matrix element f 2 , ([I]), also plays a role in the so-called "intrinsic" charm 
contribution to the nucleon spin 0. Making use of an expansion in inverse powers of the 
charm quark mass, l/m c , Franz et al. 0, correcting an earlier results [[], PJ, related the 
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charm quark contribution to the axial current to the matrix element of the flavor-singlet 
twist -4 operator, ([!]). Model estimates of flavor-singlet / 2 can thus be used to estimate 
the intrinsic charm contribution to the nucleon spin. Knowledge of this contribution 
is a prerequisite for attempts to measure the gluon polarization through open charm 
production JTT|. 

Here we report about an estimate of the spin-dependent twist-3 and 4 nucleon matrix 
elements, (|ip and (^), within an instanton-based description of the QCD vacuum. In this 
approach the QCD ground state is approximately described as a "medium" of topologi- 
cal vacuum fluctuations — instantons and antiinstantons, with perturbative fluctuations 
about them |12|, [L3l KM. This picture explains the dynamical breaking of chiral symmetry 



15] , which happens due to the fermionic zero modes associated with the individual (anti- 



instantons, and manifests itself in the appearance of a dynamical mass of the quarks, ac- 
companied by a coupling to Goldstone bosons — the pions. Building on that, this picture 
explains a host of phenomenological data on "vacuum structure" and hadronic correlation 
functions ||14|| . It also gives rise to a successful description of the nucleon as a chiral soliton, 

oo 
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in the formal limit of a large number of colors, N c 

A crucial element in this approach is the fact that the instanton medium is dilute, 
i.e. the average size of the (anti-) instantons, p ~ (600 MeV)^ 1 , is small compared to 
their average distance, R ~ (200 MeV) -1 , with p/R w 1/3. The existence of this small 
parameter makes possible a systematic analysis of non-perturbative effects generated by 
the instantons. Matrix elements of quark-gluon operators of twist > 2 have been studied 
in this approach in Ref.[18|. In particular, it was shown there that the matrix element 
of the twist-3 operator, d 2 , is suppressed relative to the twist-4 one, / 2 , by a factor of 
(p/R) 4 , and thus 

4 « /, 



twist-3 



twist-4 



In Ref.||18|| the value of d 2 was estimated to be of the order of 10 ; this prediction is 
confirmed by the recent results of the E155x experiment || . Thus, in the instanton vacuum 
the parametrically leading higher-twist effects are due to the twist-4 matrix element f 2 ■ In 
this letter we present a quantitative estimate of f 2 in the instanton vacuum. In addition to 
the flavor-nonsinglet nucleon matrix element, which appears in leading order of the 1/N C 
expansion and has been calculated in Ref.|18], we calculate here also the flavor-singlet 



one. This allows us to make predictions for the power corrections to both the proton and 
neutron spin structure functions G\, as well as to estimate the intrinsic charm contribution 
to the nucleon spin. For a detailed description of the foundations of the approach employed 



in the present investigation here we refer to Refs.[|T^, |1| 



2 Twist— 4 matrix elements from the instanton vac- 
uum 



Chiral symmetry breaking by instantons. Following Diakonov and Petrov [15, 20], the 
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dynamical breaking of chiral symmetry in the instanton vacuum can be studied in the 
large- N c limit, where it manifests itself in the appearance of a dynamical quark mass, 
M^O, accompanied by a coupling of the quarks to a pion (Goldstone boson) field, 7r a (x), 
(a = 1, . . .Nj — 1) which arises from the "bosonization" of the 't Hooft many-fermionic 
interaction induced by the instantons. The resulting low-energy dynamics is summarized 
by the effective Lagrangian 

L cS = tfj(x) [ijd- MF(d 2 )U' r5 (x)F(d 2 )] t/j(x), (6) 

where 

U^(x) = ^-^U{x)+— —U\x), U(x) = e i7Ta{x)r \ (7) 

is a unitary matrix variable containing the pion field, and F(p 2 ) are form factors, propor- 
tional to the Fourier transform of the zero mode of the (anti-) instantons, which drop to 
zero for spacelike momenta of the order of — p 2 ~ p~ 2 . It is crucial that the dynamical 
quark mass, M, is parametrically small compared to the ultraviolet cutoff, p _1 : 



Mp 



In particular, this implies that for spacelike quark momenta of the order M <C ~ 
600 MeV one may neglect the form factors and consider the quark-pion coupling as effec- 
tively pointlike. 

Nucleon as chiral soliton. The effective Lagrangian @ serves as a basis for the calcula- 
tion of correlation functions in the instanton vacuum within the l/iV c -expansion j[5|, |20fl . 
In particular, it gives rise to a picture of the nucleon as a chiral soliton [Tj|. Baryon cor- 
relation functions in the large-iV c limit are characterized by a classical pion field, which 
in the baryon rest frame is of "hedgehog" form 

UJx) = e inVap(r) , r = IxL n = -. (9) 

r 

Baryon states of definite spin/isospin and momentum quantum numbers emerge from 
quantizing the collective rotations and translations of the classical soliton 

U(x,t) = R(t)U c (x-X(t))R\t). (10) 

Since the soliton moment of inertia and mass are of order N c , the angular velocity of the 
collective rotation, 

n = -iR ] R, (11) 

as well as the linear velocity of the translation, are of order 1/N C , the quantization of the 
collective motion can be performed within a. l/N c expansion []16| . The (^-contribution 
to the baryon mass gives rise to the N—A mass splitting, which is of order 1/N C . This 
approach allows also to compute matrix elements of operators between baryon states. 
Technically, this is done by representing the matrix element as a functional integral over 
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Q— + -0 — *- 



Figure 1: Schematic illustration of the contributions to the matrix element of the quark- 
gluon operator, (|l|), in the instanton vacuum. The cross denotes the quark-gluon operator, 
the solid line the quark propagator, the dashed line the gauge field of the (anti-) instanton, 
and the circle the chirality-flipping 't Hooft interaction of the instanton with the quark 
mediated by the zero modes. The picture shows the situation for one light quark flavor; for 
Nf flavors the 't Hooft interaction is a 2Nf- fermionic vertex, which at large N c is equivalent 
to a coupling of the quark to the pion field ( "bosonization" ) . Using the equations of 
motion of the effective low-energy theory, the instanton-induced vertex shown here can 



be converted to the chirally even effective operator (|12|). For details, see Refs. |18|, 21 . 



the collective translations/rotations, which is then evaluated by expanding the integrand 
in fl. For details we refer to Refs.[|TB], |T7| . 

Twist-4 quark-gluon operator in the instanton vacuum. The instanton vacuum, being a 
microscopic model of the non-perturbative fluctuations of the gluon field, allows to estimate 
hadronic matrix elements of quark-gluon operators, such as the twist -4 operator of Eq. ([![) 
T9| , 18 ]. It is understood that the QCD operators are normalized at the scale defined by 
the inverse average instanton size, \x = p~ l . The main idea of this approach is that because 
of the diluteness of the instanton medium the gauge field operator in (Jf]) can effectively 
be replaced by the field of a single (anti-) instanton, which then interacts with the quark 
fields through its zero modes (this is schematically illustrated in Fig. p. As a result, the 
original QCD quark-gluon operator is replaced by a chirality-flipping "effective operator", 
which is to be evaluated in the effective theory where chiral symmetry is spontaneously 
broken. As has been discussed in Refs. [|T8| , |21~H , this effective operator can be brought to a 
chirality-conserving form by making use of the equations of motion of the effective theory. 
For the QCD operator of interest, tpfjaFpaipf, the resulting operator is of the form 



(12) 



It is understood here that the momenta of the quark fields are of the order p ~ M <C p _1 . 
Here Z^"^ denotes the matrix element of the original chirality-flipping effective operator 
in a quark state, given by 

;r rk = h(- P 2 )\^=M*, (is) 

where p 2 ) denotes the quark loop integral defined by (p and k are Euclidean momenta, 
with p 2 = —p 2 ) 



Hp 2 



al 4 k F{p-k)g{k) 
(2tt) 4 (p-kf 



■!§<*><«.. 



cos( 



(14) 
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in which are the Gegenbauer polynomials (4-dimensional spherical harmonics) of 
argument cosO = (k ■ p)/(\k\\p\), and the function Qik) describes the Fourier transform 
of the (dual) field strength of a single (anti-) instanton in singular gauge (see Ref.jnj for 
details): 



g(k) 



32vr z 



t 1 



t 



\k\p, (15) 



with K n (t) the modified Bessel functions of the second kind. The numerical value of /2 Uark 
in the limit Mp — > is ||T8 

3 



quark 



0.6. 



(16) 



The instanton-induced effective operator ([12]) has the form of the axial current oper- 
ator, but with an additional contracted derivative acting on the quark fields. Note that 
the coefficient, /2 Uark ; is °f order unity in the packing fraction of the instanton medium, 
p/R. Technically speaking, this happens because the quark loop integral 1%, (f4]), contains 
a would-be quadratic divergence, which is regularized by the instanton form factors and 
cancels the overall factor p 2 in front of the integral. This implies that the hadronic matrix 
elements of this operator will also be parametrically large in the instanton packing frac- 
tion. We remark that for the twist-3 operator in the matrix element, d 2 , (0), the effective 
operator would be proportional to (Mp) 2 ~ (p/R) A and thus parametrically suppressed. 

Nucleon matrix elements: 1/N C expansion. The nucleon matrix elements of the effective 
operator, fll2|) , can be computed using the same techniques as have been used in the 
calculation of nucleon matrix elements of vector and axial vector current operators (for a 
review, see |JL7|| ). In particular, we shall exploit the close relation of the effective operator, 
(P^), with the axial current operator. It is known that within the chiral soliton picture of 
the nucleon the isovector and isoscalar component of the nucleon axial coupling constant 
appear in different orders in the expansion in powers of the angular velocity of the collective 
rotation of the soliton, Q ~ 1/N C , fllTI). We therefore need to consider separately the 
isovector and isoscalar components of the twist-4 nucleon matrix element f 2 , (P]), 



f (3) 
J2 



J2,u ~ J2,d, 



AO) 

J2 



J2,u + f2,d 



(17) 



(the matrix elements on the R.H.S. refer to the proton). We restrict ourselves for the 
moment to the SU(2) flavor group (u,d), the extension to SU(3) will be discussed below. 
For the isospin components defined by (|I7]) the ^-counting can immediately be inferred 
from that of the corresponding axial coupling constants of the nucleon: 



M 2 N f 2 



Nr. 



1. 



The isovector component of the spin-dependent twist-4 matrix element is leading in the 
1/N C expansion. 

(3) 

Isovector matrix element / 2 . We now make a quantitative estimate of the isovector 

(3) 

twist-4 nucleon matrix element, / 2 • This can be done by considering the hypothetical 
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limit of large soliton size, in which one can perform an expansion of the matrix element 
in gradients of the classical pion field, @. The result for in this approximation can 
immediately be written down by analogy with the well-known result for the isovector axial 



coupling, gf ] |6] 



(3) 
9 A 



AN C M 2 J 
9 



d x tr 



-iT a Ui(-K)d a U c (-x) 



(19) 



where J denotes the (Euclidean) quark loop integral 

d ±k F\k) \F\k) - k 2 F{k)F'(k) 



J 



(2tt)< 



k 2 + M 2 F 4 (k) 



(20) 



(Note that the PCAC relation relates this integral to the pion decay constant, 4N C M 2 J 
F 2 .) The corresponding result for the isovector matrix element of the operator (|T2"D is 



4AT c M 2 / 2 quark J 1 
9 



d x tr 



-IT 



a [/t(x)9 a C/ c (x) 



(21) 



where 



Ji 



which differs from J, 



d*k k 2 F 2 (k) F 2 (k)-k 2 F(k)F'(k) 



(2tt)< 



k 2 + M 2 F\k) 



[22] 



by an additional power of the quark virtuality, k 2 , in the 



integrand. For the ratio of / 2 to g^' we thus obtain 



,(3) 



(3) 
9 A 



f: 



quark j 



J 



(23) 



With the standard parameters of the instanton vacuum, p = (600 MeV) 1 and M 
350 MeV, we find 



/.quark j 

h •h 

J 



0.49 p~ 2 . 



With the experimental values M N = 940 MeV and g A 1 = 1.25, Eq.(^) thus gives^ 



(3) 



-0.25 



(24) 



(25) 



The integral J\ in the gradient expansion of / 2 , (^2|), contains a would-be quadratic 
divergence, i.e., it is parametrically of the order Ji ~ p~ 2 , contrary to the integral J in 
g^\ which is depends only logarithmically on p~ x . Thus, in J\ the dominant contribu- 
tion comes from quark momenta of the order of the UV cutoff, k ~ p~ x . In principle 

1 In the numerical estimate quoted in Ref.jl8]], = —0.11, some contributions to the gradient expan- 
sion of the matrix element were missed. Hence the differences in the numerical values. 
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the representation of the instanton-induced effective operator by the local operator, (12), 
is accurate only for quark momenta of the order k ~ M <C p _1 ; for momenta of the 
order JT X one should take into account the momentum-dependence (non-locality) of the 
instanton-induced effective operator. One should keep in mind that, anyway, the precise 
form of the UV cutoff of the effective low-energy theory depends on the details of the 
approximations made in the description of the instanton medium; only the gross features 
can be assumed to be generic. Thus, results for quantities given by "quadratically di- 
vergent" integrals, such as , should generally be regarded as rough estimates (±50%). 
Whether or not to include the form factors in the effective operator is just part of this 
larger uncertainty. We have verified that using instead of (O) the "exact" non-local op- 
erator —^7/375 [h(d 2 ) / F (d 2 )] d 2 ijj changes the numerical result for / 2 by less than 10%, 
so using the local approximation to the operator seems completely justified. 

Isoscalar matrix element f^. The calculation of the 1 /^-suppressed isoscalar com- 
ponent of the twist— 4 matrix element, / 2 , is somewhat more involved, since it requires to 
take into account the time dependence of the saddle-point pion field, (|10"1), to first order 
in the angular velocity of the soliton, Q, QTTJ) . An important question is the degree of 
the would-be ultraviolet divergence of the isoscalar nucleon matrix element, f%. This 
question can be investigated by gradient expansion. Expanding the average over quark 
fields of the flavor-singlet version of the operator (|12j), 



& 1 

L Flavor 7/375^ 



ijd - MF(d 2 )U^F(d 2 ) 



(26) 



in gradients of the (space- and time-dependent) background pion field, one finds that the 
functional trace is at most logarithmically divergent. This implies that the matrix element, 
which is obtained by substituting in fl2"6f) the slowly rotating pion field (0), and expand- 
ing to first order in the angular velocity, Q, can also be at most logarithmically divergent. 
This circumstance is very fortunate, as it allows us to neglect the instanton-induced form 
factors in the quark-pion coupling in the Lagrangian, (P), and simply apply an external 
UV regularization in the form of a Pauli-Villars subtraction. (If the matrix element were 

(3) 

quadratically divergent, as / 2 IS > the result would have been strongly dependent on the 
specific form of ultraviolet cutoff applied.) With this simplification the calculation of the 
matrix element of the instanton-induced effective operator (0) becomes completely anal- 
ogous to that of the isovector axial coupling constant, [H, ID, anc ^ can be performed 
using the same techniques. 

A new feature in the calculation of the matrix element of the instanton-induced oper- 
ator ( |T2"D compared to that of the axial current, ^7^75^, is that with (|I2|) contributions 
of order Q l can appear from the action of the derivatives contained in the operator on 
the time-dependent isospin rotation matrices, R(t). Such contributions were discussed in 
Ref. [2_3]| in the context of the calculation of the isovector unpolarized quark distribution in 



the nucleon, and we refer to this article (in particular Section 3.1) for a detailed description 
of the Q- (1/N C -) expansion. After performing the fi-expansion one expresses the remain- 
ing Green functions in the static background pion field, (|]), in a basis of eigenfunctions of 
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the quark single-particle Hamiltonian in the: 



H = -i 7 VV< + 7°Mf/ 75 , (27) 
H\n) = E n \n). (28) 



The result for the isoscalar matrix element then reads 

f (o) Njr rk J i V V 1 



■\l-s 1/ n m 

E n — E 
in 



non— occup. occup. 

(n|<T 3 (E^ - p 2 )|m) (m\r 3 \n) + (n\r 3 \m) {m\a 3 {El - p 2 )\n) 



-777 J2 En(n\a 3 T 3 \n) 



21 

occup. 



(29) 



Here / is the moment of inertia of the classical soliton, see Ref. ||16||. The first term in 
the braces is a double sum over quark levels; it coincides with the expression for g$ up 
to the insertions of the operator E 2 — p 2 in the single-particle matrix elements, which 
come from the contracted derivative in the operator (|i~2"|). The second term represents the 
abovementioned "new" contribution which results from the action of the derivative in the 
operator on the rotational matrices in the collective quantization. Up to the factor E n 
inside the single sum over quark levels this contribution is identical to the expression for 

(3) 

the isovector axial coupling constant, g A . 



In Eq. (p9|) the first term (the double sum over quark levels) turns out to be UV-finite 
and does not require regularization. The second term (the simple sum over levels, related 
to the isovector axial coupling constant), as it stands, contains a linear divergence when 
regularized e.g. with an energy cutoff. At the same time one observes an anomaly-type 
phenomenon in the sense that the sum over occupied quark levels is not equal to minus 
the sum over non-occupied ones, as it should be on grounds of the analyticity proper- 
ties of the single-particle Green functions. However, both the linear divergence and the 
"anomaly" disappear when performing a Pauli-Villars (PV) subtraction, leaving behind 
a finite sum with usual behavior when summing over non-occupied instead of occupied 
states, in agreement with the above statement that the matrix element is at most loga- 
rithmically divergent. A similar phenomenon was encountered in the calculation of matrix 
elements of twist-2 operators of spin > 2 in the chiral quark-soliton model in Ref. |25| . 



In that case, too, the superficial power divergences and the "anomalies" simultaneously 
disappeared after PV subtraction. Note that the PV subtraction, contrary to regulariza- 
tion methods based on an energy cutoff, does not spoil the analyticity properties of the 
regularized sums. 

The remaining PV regularized sums can be performed numerically, using the Kahana- 
Ripka method of diagonalizing the single-particle Hamiltonian in a basis of free quark 
states. The value of the PV cutoff is determined by fitting the pion decay constant f25fl , 



9 



Mpv = 557 MeV. For the numerical estimate we use the variational soliton profile of 
Refs. [|16| , 25f1 . We obtain a numerical value of the flavor-singlet twist -4 matrix element of 

/ 2 (0) = 0.01. (30) 
The flavor-singlet matrix element is more than an order of magnitude smaller than the 



flavor-nonsinglet one, (25]), which reflects the fact that the former is given by a "logarith- 
mically divergent", the latter by a "quadratically divergent" expression. 

Comparison with results of other approaches. The instanton result for the isovector 

(3) 

twist -4 matrix element, / 2 ■> agrees both in sign and in magnitude with the QCD sum 
rule results of Refs. |2T)|, [2"Tf . However, we disagree with these authors on the sign of the 
isoscalar matrix element (for a critical discussion of the QCD sum rule calculations, see 
Refs.[^]). We also disagree in the sign of the isovector matrix element with the bag model 
estimate of Ref. 0; however, this model can hardly claim to give a realistic description of 
the quark-gluon correlations giving rise to the twist-4 matrix elements; e.g, it does not 
respect the QCD equations of motion. 

Generalization to the SU(3) flavor group. So far we have calculated the twist-4 matrix 
elements assuming the SU(2) flavor group (u, <i-quarks only). The results can easily be 
generalized to the SU(3) case (including also s quarks). Here we consider the simplest 
scenario of perfect SU(3) symmetry, i.e., zero strange quark mass. In the chiral quark- 
soliton model the leading-iV c calculation gives a simple relation between the SU(3) triplet 
and octet couplings, as the differences between them appear only due to the different 
states of collective rotations of the soliton. The ratio of the triplet and the octet part of 
the nucelon matrix element is therefore given by the ratio of the matrix elements of the 
respective Wigner D-functions in the rotational state corresponding to the nucleon spin, 
isospin and the hypercharge, |rot) = \S = T = 1/2, S3 = 1/2, T 3 = 1/2) (for details see 

(31) 



if (rot I £>» 3 I rot) 3 



( rot 


D 8 

-^83 


rot ) 


(rot|D§3 


rot) 



/( 3 ) ' (rot I rot) " 7' 

From these relations we obtain the following numerical results for the twist-4 matrix 
elements in the SU(3) symmetric case: 

(32) 
(33) 
(34) 



AO) 

J2 


SU(3) — 


0.01 


A3) 
J2 


SU(3) — 


-0.25 


A8) 
J2 


SU(3) = 


-0.11 



3 Power corrections to proton and neutron spin struc- 
ture functions 

With the results obtained in the previous section we can now estimate the 1 / Q 2 -corrections 
to the first moments of the proton and neutron spin structure functions, Gi'°. The expres- 
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sions obtained from the operator product expansion of QCD at tree level to order 1/Q 2 
are HQ 

* dxG^ix, Q 2 ) = l -a^ n) + ^ (a p 2 (n) + Adf n) + 4/f W ) . (35) 

Inclusion of radiative corrections would lead to a logarithmic Q 2 -dependence of the coef- 
ficients. Here denote the spin-dependent matrix elements of the twist-2 spin-(n + 1) 
axial vector operator; the proton and neutron matrix elements (which include the quark 
charges) are given in terms of the SU(3) singlet, triplet and octet components as 

^ (n) =4 a(3) +^ a(8) +^ (o) ' ( 3e ) 

with +(— ) for proton (neutron). We quote the expressions for three light quark flavors 
(u, d, s). For n = the a n coincide with the axial coupling constants of the nucleon: 

4 3) = 9 ( 1\ 4 8) = gf, 4 0) = 9 f- 07) 

Of the power corrections the term proportional to are the target mass corrections, 
while the terms proportional to and Jf represent the dynamical higher-twist cor- 
rections; the proton and neutron matrix elements are defined in analogy to Eq. (|36"D . 

We now evaluate the dynamical power corrections using the instanton vacuum results. 
Following the basic philosophy of the instanton vacuum we put the parametrically sup- 
pressed twist-3 spin-2 matrix elements to zero (these matrix elements were estimated in 
Ref . |18|| to be of the order of ~ 1% of the twist-4 ones). For the twist-4 matrix elements 



ff we use the results (|3~2"D-(|3"4T), which imply 

fP = -0.046, ft = 0.038. (3* 



To evaluate the twist 2 contribution we use the GRSV 2000 parameterization |29| ( "stan- 
dard scenario"), which includes the radiative corrections to the coefficient functions in 
NLO.0 The results are shown in Fig.^. One sees that the dynamical twist-4 contribution 
has a rather small effect on the structure functions down to Q 2 of the order of 1 GeV 2 , 
in particular in the neutron. Note that the quantitative details may change when SU(3) 
symmetry breaking effects are included in the twist-4 contribution. 



2 Note that there is a mistake in Ref.|l| concerning the coefficients in this formula, as was noted in Q; 
see also Ref. [|J . 

3 To evaluate the target mass corrections we use instead of the twist-2 matrix element a^ n ^ the third 
moment of the twist-2 part of G^ n ^ calculated from the GRSV 2000 parametrization; the difference 
between the two is irrelevant in the present context. 



11 



Proton 




Neutron 
-0.04 i — ■ — ■ — ■ — ■ — ■ — ■ — ■ — - 




1 5 g 2 /GeV 2 10 



Figure 2: The Q 2 -dependence of the first moment of the polarized structure function G\ 
for the proton (top) and neutron (bottom). Dashed lines: Twist-2 contribution according 
to the GRSV 2000 NLO parameterization (standard scenario) p§| . Dotted lines: Sum 
of twist-2 contribution and target mass corrections. Solid lines: Total result, including 
also the twist-4 contribution due to the matrix elements f2, as estimated in the instanton 
vacuum. (In the case of the neutron the target mass corrections to the first moment are 
very small, so we show only the pure twist-2 contribution.) 
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In Fig.^| we show the result for the first moment of the difference of proton and neutron 
structure functions, G\ — G" (the Bjorken sum). This quantity is of particular interest, 
since the Bjorken sum rule is a rigorous prediction of QCD, and the radiative corrections 
have been calculated to order [[Kj. We again use the GRSV2000 NLO parametrization 
to estimate the twist-2 part and target mass corrections. The power corrections now 
receive contributions only from the flavor-triplet twist-4 matrix element, / 2 ■ Also shown 
in Fig.^| are the experimental results obtained in the analyses of the SMC E143 [32], 
and E155 |33| experiments. One observes that the relatively large twist-4 correction 



obtained from the instanton vacuum improves the agreement of the theoretical prediction 
with the data; however, the present experimental errors are too large to allow for definite 
conclusions. 



0.2 



r 



0.15 



Bjorken SR 



p-n < 




5 gVGeV 2 10 



Figure 3: Same as Fig.|2|, but for the difference of proton and neutron structure functions 
(Bjorken sum rule). Also shown are the experimental results quoted by the E143 collabo- 
ration at Q 2 = 2, 3 and 5 GeV 2 [§2J, and by the SMC |3J and E155 
Q 2 = 5 GeV 2 . (The abscissae of the data points at Q 2 
in order to separate them in the plot.) 



collaborations at 
5 GeV^ have been shifted slightly 



4 Intrinsic charm contribution to the proton spin 

The isoscalar component of the twist-4 matrix element, f^, (|T|), plays a role in the nucleon 
spin structure functions aside from power (l/Q 2 -) corrections, namely in determining the 
intrinsic charm contribution to the nucleon spin. Up to power corrections, the first moment 
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of the proton spin structure function for three light flavors (u, d, s) is given by [cf. Eqs. (|35| 
and ©] 



When the charmed quark is included, its contribution can be expressed as a modification 
of the SU(3) flavor-singlet axial coupling constant in (]39|): 



gf = g%w* + 2g { 2\ (40) 

where g^ught contains the S'C/(3)-singlet contribution from the light flavors only, and 
is defined as 

(p,s|cy/j75cb,s> = -2gys p . (41) 

In the proton, where the charmed quarks occur in the form of virtual cc pairs, the charm 
quark axial current operator can be approximated by an operator involving only light 
quarks (u,d,s) and the gluon fields, using the heavy-quark mass (l/m c -) expansion to 
integrate out the charm degrees of freedom. Refs.|| [7], correcting earlier results of Refs.|| 
10| , quote the following result for the nucleon matrix element of the charm axial current 



to order l/m 2 c : 



(P, S|c7/37 5 c|p,s) = ' 2 (p,s\ 22 ^na F l3a1pf\p,s), (42) 

izirm c f= u ,d,s 



from which follows that 

127r m 2 2 



gT = ~-^fr- (43) 



c 



Using our result for the SU(3) flavor-singlet matrix element, /2°\ Eq. (|3"2]) , we can now esti- 
mate the charm axial constant, g^. With a s (fi = m c ) = 0.39 and m c = 1.15 . . . 1.35 GeV, 
we find 

g^ ~ -1 x 10" 4 . (44) 

The charm quark contribution to the proton spin is very small, mostly because of the 

■(o 

2 



smallness of the flavor-singlet twist-4 matrix element, //° 



5 Summary 

The instanton vacuum with its inherent small parameter — the packing fraction of the 
instantons, p/R — implies a parametrical (and numerical) hierarchy of the spin-dependent 
twist-3 and 4 matrix elements: twist-3 twist-4. The unusually small value for the 
twist -3 matrix element ~ 10 -3 obtained from the instanton vacuum [18] agrees well 
with the precise measurements of G<i in the E155x 0] experiment (and, also, the older 



E143 [p2j and E155 |33| results). Note also that recently revised lattice calculations of c?2 



34], properly accounting for operator mixing effects, confirmed the small value predicted 



by the instanton vacuum. 
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In this paper we have presented numerical estimates of the twist -4 matrix elements 
/2, including also its flavor-singlet component. Our approach predicts a big numerical 
difference between the flavor-nonsinglet and singlet matrix elements, {f^l 3> If^l, in 
contrast to QCD sum rules, which tend to give values of the same order of magnitude. 
Phenomenologically, this would imply much larger power corrections to the Bjorken than 
to the Ellis-Jaffe sum rule, a prediction which can in principle be tested experimentally. 
We have also pointed out that the small value of f^ implies a very small intrinsic charm 
contribution to the nucleon spin. 

The approach laid out here can be extended in many ways; for example, one can 
compute also the matrix elements determining power corrections to higher moments of the 
structure functions and restore the the ^-dependence of the twist-4 contribution. This 
would allow a more direct comparison of the higher-twist corrections with experimental 
data. 
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